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We consider the effect of new physics on the branching ratio of Bg l^l^^ where 
I = e, fi. If the new physics is of the form scalar/pseudoscalar, then it makes no con- 
tribution to Bg l^l^^, unhke in the case of Bg , where it can potentially 
make a very large contribution. If the new physics is in the form of vector/axial- vector 
operators, then present data on i? — > [K, K*)l^l^ , does not allow a large enhance- 
ment for B{Bs l^l'^). If the new physics is in the form of tensor/pseudotensor 
operators, then the data on B ^ {K, K*)l^l~ gives no useful constraint but the 
data on i3 ^ K*'y does. Here again, a large enhancement of B{Bs — > much 
beyond the Standard Model expectation, is not possible. Hence, we conclude that 
the present data on 6 — > s transitions allow a large boost in B{Bs l^l^) but not 
in B{Bs 1^1-7). 

The quark level transition b sl^l~ can lead to a number of important flavour changing 
neutral current (FCNC) decays in B mesons. Among them are the semi-leptonic modes B 
{K, K*)l'^l~ , the purely leptonic mode Bg l^l^ and the leptonic radiative mode Bg — >■ 
The relationship between the semi-leptonic and purely leptonic modes was discussed 
in ll||. It was shown that, if the new physics occurs in the form of vector /axial- vector 
operators, then present data on semi-leptonic branching ratios [2,131 constrain the branching 
ratio for Bs to be of the same order of magnitude as that of the Standard Model 

(SM). On the other hand, if the new physics operators are in the form of scalar/pseudoscalar, 
then the semi-leptonic branching ratios do not lead to any useful constraint on the rate for 
the purely leptonic mode. Hence a large enhancement of Bg l^l^ is possible only if the 
new physics is in the form of scalar/psuedoscalar operators. Tensor/pseudotensor operators 
do not contribute to Bg . In this letter, we examine the relation between the rates 

for effective 6 — s transitions and Bg — > l^l'j. 

In the SM, the decay Bg — > has small branching ratio due to helicity suppression. 
The radiative decay Bg l^l^^ is free from helicity suppression due to emission of a photon 
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in addition to the lepton pair. Thus the branching ratio for this leptonic radiative mode 
is much higher than that for the purely leptonic mode despite an additional factor of a. 
Because of this higher rate, this mode will be an important probe of 6 ^ sl^l" transitions 
which will be studied at present and future experiments. The decays Bg —>■ l^l^j have been 
studied in several papers 0,0,0,0, 2,0, [3 within the framework of SM. The effective new 
physics Lagrangian for b sl^l~ transition is the sum of three terms: vector/axialvector, 
scalar/pseudoscalar and tensor/pseudotensor. The first two terms can arise both via penguin 
and box diagrams but the last term arises only via the penguin diagram for 6 ^ S7, in which 
the real photon is replaced by a virtual photon coupling to a lepton-antilepton pair. In [^,0|, 
the effective h sl^l^ interaction was dressed with an on-shell photon in all possible ways. 
Helicity suppression is operative for the case where the photon is emitted from the final 
lepton and the resultant amplitude is proportional to the lepton mass and is negligible. For 
the case where the photon is emitted from the internal lines of the b ^ s loop transition, the 
amplitude is suppressed by factors ml/m^ and is also negligible. The main contribution to 
the Bg amplitude comes from the diagrams where the final state photon is emitted 

from either 6 or s quark in the effective b — sl^l~ interaction. With this procedure, the SM 
prediction for B [Bg e^e^'-f) is calculated, in to be about (2 — 7) x 10~^, with the 

rate for Bg fi^fi~'y being a little lower. 

In ref.jll a higher value of branching ratio for Bg — > /+/~7 is predicted within SM. This 
higher value is due to a different parametrization of the form factors /y, /a, frv and frv- 
The parametrization of form factors in is based on QCD sum rules whereas in ^1 it is 
based on light front models. But ref.j3| uses the parametrization based on perturbative QCD 
methods combined with heavy quark effective theory In ref.jl^, it was argued that 

there are additional contributions to the Bg — > l~^l~'y amplitude. The most important one 
comes from the case where the real photon is emitted from the b s loop transition and the 
virtual photon, which pair produces the leptons, is emitted from the initial quarks. Due to 
this additional amplitude, the SM prediction for B {Bg — > e"*"e~7) in jl^, is about 2 x 10^*^, 
with the branching ratio for Bg — > /i+/i^7, being a little smaller compared to Bg e"'"e~7. 

In the present calculation, we are interested on how the current data on b s transi- 
tions, due to the effective interactions b sl^l~ and b 57, constrain the new physics 
contribution to the leptonic radiative decays Bg l^l^^y. 

As mentioned earlier, new physics in the form of scalar/pseudoscalar operators can give a 
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large enhancement to the leptonic decay mode Bg . The question then follows: What 

is the effect of these operators on the leptonic radiative modes Bs Unfortunately, 
scalar/pseudoscalar operators do not contribute to Bs I'^l^-y. The photon has J = 1. 
Hence the pair also must be in J = 1 state so that the angular momentum of the 
final state can be zero. However, by Wigner-Eckert theorem, the matrix element {l~^l~{J = 
^)\K9s + fl'p75)^|0) is zero. This result also follows from direct calculation, as we illustrate 
below. 

We parametrize the scalar/pseudoscalar operator for b — > sl^l~ transition as 



Gp j 



a 



Lsp{b^sl+l ) = -j=\^-^-^] s{gs + gp-f5)bl{gs + 9p75)1- (1) 



The matrix element for Bs —>■ l^l 7 is given by 
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M{B ^ /+/-7) = -j= ) [gs{^ \sb\ Bsip)) + gp{i |s75&| Bs{p))] u{Pi)i9s + 9p1>{Pi) 
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To calculate the matrix elements of the quark operators in the above eciuation, we need to 
first consider the following vector and axial-vector matrix elements 

{^{k)\-s^^b\Bs{p)) = ee,,p^e*^pfk'^Mq')/mp^, 
(lik) |S7^75&| Bs{p)) = -le [t*^{p ■ k) - (e* ■ p)k^] fA{q')/mB^, (3) 

where q = Pi+Pj- Dotting the above equations with the momentum of Bs meson p^, we get 
the scalar and pseudoscalar matrix elements to be identically zero, 

{^(k) \sb\ Bsip)) = = (lik) \s^,b\ Bsip)). (4) 

That this amplitude vanishes, was also demonstrated in So, even if a large enhancement 
of Bg l^l^ is observed at LHC-b 3| due to new physics operators in scalar/psuedoscalar 
form, there will be no corresponding enhancement of Bs l^l^^y. 

A legitimate question to ask at this stage is: Is it possible to have a large enhancement 
of Bs l^l^'-y for any type of new physics operator? Here we consider vector/axial- vector 
operators and tensor/pseudo-tensor operators one at a time and examine their contribution 
to Bs given the current experimental results on the b ^ s transitions. 

First we will assume that the new physics Lagrangian contains only vector and axial- 
vector couplings. We parametrize it as 

LvAib stl~) = ^ ^(9v + 9Alb)liJ,b lig'y + g'Al5)Yl, (5) 
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where g and g are effective couplings wliicli cfiaracterise tfie new pliysics_^ 

Tlie vector and axial-vector matrix elements are shown in Eq. Q 
dence of the formfactors is parametrized as P. [lfll|. 

fBs^Bs 



where i = V,A, TA, TV and the parameters /? and A are given in Table I. 



Table I: Parameters for the form factors 



Parameter 


fv 


frv 


fA 


Ita 


P{GeV-^) 


0.28 


0.30 


0.26 


0.33 


A{GeV) 


0.04 


0.04 


0.30 


0.30 



The depen- 
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The calculation of decay rate gives 



30727r4s,^ 



where (z = V,A) are the integrals over the dilepton invariant mass {z = q'^/m%^ 
are given by 



dz ■ 



z{l-zY 



(7) 
They 

/o [(A,/mBj + 0.5(1 -z)f 
Here we have neglected the lepton masses in comparison to as we are only considering 

/ = e, /i. We will work under this approximation throughout the paper. 

In order to put bounds on B^p [Bg we need to know the values oi gy (^g'y + g'J^ 

and g\ (^g'y + g'£j . For this we will have to consider the semi-leptonic decay modes B — »• 

{K, K*)l^l~ . The values of these quantities were calculated in (ij, 



9l{9^ + 9l) = (1.36l°iDxlO 
9\{9v + 9l) = {G.rGttit) X 10 



2 
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These values were calculated under the assumption that B^p [B — > {K, K*)l^l^] = 
Bexp [B {K, K*)l'^l~] i.e. the experimentally measured semi-leptonic branching ratios 
are saturated by the new physics couplings. Putting these values in Eq. (0), we get 



Bnp (b, i^r-f) = 2m1:^ofe x lO" 



(10) 
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Therefore the upper bounds on the branching ratios are, 

Bnp [Bs l^l'l) < 2.90 X 10"^ at la 

Bnp {Bs l^l'l) < 4.58 X 10"^ at 3a. (11) 

These values are of the same order of magnitude as SM prediction. Thus we see that 
we can't boost Bjyp (Bg l^l^l) above its SM prediction even after assuming that the 
contribution to the decay rate is totally due to new physics. The fact, that the experimentally 
measured values of the semileptonic branching ratios B{B {K, K*)l^l~) are close to their 
SM predictions, doesn't allow Bj^p [Bg l~^l~'y) to have a value much different from its SM 
predictions if the new physics responsible for this decay is of the form vector/axial- vector. 
A more stringent upper bound is obtained if we equate the new physics branching ratio to 
be the difference between the experimental value and the SM prediction. In fact, this upper 
bound is consistent with zero at 1 a and is Bjyp^Bg I'^l^l) < 2 x 10~^ at 3 a. Thus we 
can't boost B^p {Bs — > l^l^'j) much above its SM prediction if new physics is of the form 
vector/axial- vector. 

We now consider new physics interaction in the form of tensors/pseudo-tensor operators. 
We parametrize this Lagrangaian as, 

L^^b sin-) = ^ (^^^ (^^^ rsa,.q-'{9TV + 9TAl,)h l^l- (12) 

The necessary matrix element for Bs —>■ /^/~7 is given by, 

{l{k)\sta,,q''b\Bs{p)) = -ee^,p.e*VA:VTv(g'), 
{l{k) I w^.75g"&| Bs{p)) = -le [e;(p • k) - (e* ■ p)k,] frAiq'). (13) 

The dependence of the formfactors is given in Eq. ©. The calculation of decay rate gives, 
Tnp {Bs - /+/-7) = (^^^) [a'TvPhlTV + 9'taP'taIta] , (14) 

where /j (z = TV,TA) are the integrals over the dilepton invariant mass {z = q^/m?B^). 
They are given by 

/■I il - zY 
Ii= dz ^- — 2 . (15) 

A4mf/ms2) ;z[(Ai/mBj+0.5(l-2)]' 

Here again we have neglected the lepton masses everywhere except in the lower limit of the 
integral Jj due to the presence of term I/2; in the integrand. Thus we expect larger value for 
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Tnp {Bs — > e"'"e~7) in comparison to T^p {Bg — > due to presence of term Inz in the 

expression of decay rate. We need to know the values of g^y and g^^ in order to obtain the 
upper bound on B^p (Bg For this we will consider first the semi-leptonic decays 

B {K, K*)l+l^ and then the radiative decay B K*^. 

In order to obtain bounds on gf^v, we will have to consider the process B Kl^l~ . The 
necessary matrix element in this case is [IJ, , 

{K{pk) \ma,,q''b\ B{pb)) = -. ^ ^q\pB + Pk)^^')- W 

{nip + niK* ) 

In above equation we have dropped a term proportional to as it will give rise to a term 
proportional to (mi/m^)^ in the decay rate. The dependence of the formfactor is assumed 
to be 

MO) 



Ml') = I, ' \ 2V (17) 

(1 - q^/mj^) 



The calculation of decay rate gives, 



' G\.m%,\ ( a ^ ^ 



r„.(B^A-^r)^,|,^^j^^j /|(o)/- (18) 

where IpK is the integral over the dilepton invariant mass [z = q'^/m%J. This integral is 
given by 

where 0(^) = (^ - 1 - ey - 4p with k = niK/mB and the limits of integration for z are 
given by Zmin = 4:m'i/m% and Zmax = (1 - /c)^. 

Here again we make the approximation T^p = Tpxp ■ Under this approximation we get 
from Eq. (HH), 

2.35 [/+(0)]' ■ ^ ^ 

In order to obtain bounds on o^a, we will have to consider the process B —>■ K*l^l~ . The 

necessary matrix elements in this case are |14 Il5|. 

{K*{pu) \sia^,q''h\ B{pb)) = iT^{q^)t^,p.t*'' {ps + Pk)''{PB - Pk)" . 
{K*{pk)\sta^,q''^,b\B{pB)) = T2{q^){ml - ml,)el + T;{q^){e* ■ pb){pb + Pk*)^.. (21) 

Here again we have dropped the terms proportional to q^. The dependence of the form- 
factors is assumed to have the from, 



m) 

[I - q^/ml)' 



T-M') = n (22) 
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where i = 1, 2, 3. 



The calculation of decay rate gives, 
Tnp{B K*tl-) 



1927r3 



a 



(23) 



where I^^* {i = TV,TA) are the integrals over the dilepton invariant mass {z = q'^/m^J. 
They are given by 



tBK* 
J- TV 



tBK* 
^TA 



dz 



z{l-zf 
dz 

zHi-zy 



3/2 



Z 



{1-k 



*2\2 



8k 



1 ^{l-k*'){l-z 



20(z 
k*' 



2z + 



'1 -k*'-z 



Ak 



*2 



*2 



Ak 



*2 



(24) 



with k* 



■niB 



■. Here we assumed T2(0) ~ ^3(0). (j){z) and Zmax are the same as in the case 
of B ^ Kl^l~ with k replaced by k* . For B —>■ K*l^l~ , the experimental branching ratio is 
given with a lower cut on the di-lepton invariant masSjin^;- > 0.14 GeV, in order to supress 
background from photon conversions and 7r° e+e"7 [3]. We use this cut as the lower limit 
of integration for z. In all previous kinematic integrals, Zmin, the lower limit of intergration 
for z = q' /m\ taken to be the theoretical minimum Amf/m^. The kinematic integral, 
I^^* in Eq. (f^ , contains a 1/z' term which comes from the propagator of virtual photon 
pair producing a lepton anti-lepton pair. At very small values of z, this term dominates 
the integral and makes it very large. However, experimentally the lower limit on q' is much 
larger than the theoretical lower limit. Therefore, in calculating the bounds on new physics, 
the lower limit of in the theoretical calculation should be the same as the experimental 
lower limit. 

Under the assumption T^p = Texp and using Eq. (O!^ . we get 

Bexp{B K*l+l-) X 10^ - 1.37/#f *T2(0)^| 



2 

9ta 



Itv 



1.37/#f*^2 



(25) 



In our calulation we take the formfactors to be 



*T|(0) 

Q, /t(0) = 0.355l°:°^L TM = 0.379ir'' 



r2(0) = 0.379t o ols. The experimentally measured values of the branching ratios are 



-0.045 

Bexp{B KI+1-) = (4.8tJ:^±0.3±0.1) x 10"^ and BsxpiB ^ K*l+l 
0.2) X 10^^. Adding all errors in quadrature, we get g^y 
The best fit values for g^^ turn out to be negative and very small {O ~ 10^® 



l.Q3tlil X 10-2 for / 



11.5lii±0.8± 
= e,fx. 
The fact 
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that these come out to be negative means that the semi-leptonic decay rates can not be 
explained purely in terms of tensor/pseudo-tensor operators. Imposing the condition that 
g^Y and g^j^ should be non-negative, gives us the conditions 

g^A = 0, and g'^y = '^■Q^-oio x 10"^ for Z = e, /x. (26) 

The branching ratio for Bg due to Lt is, 

B^P [b, ^ /+/-7) = [3.15 Itv9tv + 3-81 Ita9ta] fk x 10"^ (27) 

Substituting Jb, = 240 ± 30 MeV and the values of g^y and gj,j^ in Eq. we get 

Br,p {Bs ^ e+e-7) = imt'oil X 10-' 

Bnp (Bs ^ yU+/i-7) = 6A5tlit X lO-l (28) 



Therefore the upper bounds on the branching ratios are, 



Bnp [Bs ^ e+e-7) < 2.57 X 10-^ 

Bnp {Bs /iV"7) < 8.69 X 10"^ (29) 



at Icr and 



Bnp [Bs 6+6-7] < 3.89 x 10"^ 

Bnp (Bs /iV"7) < 1-32 x 10"' (30) 



at 3(7. These branching ratios are about 40-50 times greater than the predictions in 0, 
0]. Thus the data on semileptonic decays allows an enhancement of one to two orders of 
magnitude in B^p {Bg — » /+/~7) if new physics interactions are of type tensor/pseudo-tensor. 

Here we note that b ^ S'j transition also has a tensor operator and we consider the 
constraint on the tensor/pseudotensor contribution to Bs l^l'j from the experimentally 
measured value of the branching ratio of B ^ K*'y. For this we consider the quark level 
interction 6 — > S7. We parametrize new physics effective Lagrangian for 6 — > 57 as 

L{b .7) = (li^^) ^^f^'^^'^idTV + gTAl,)he^^^^. (31) 

where e'^'^-''^ is the polarization vector of the photon and q'^ is its momentum. Replacing 
g(7)/i i-,y i^e/q^)l'^^l gives rise to 6 tensor/pseudotensor operators. Thus the present 
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experimental limit on i? — ^ ii'*7 leads to a bound on Bg /"'"/~7 arising from new physics 
operators of tensor/pseudotensor form. 
The amplitude for B K*'y is given by, 

A{B K*^) = (^^^ (l^^) ^^^^'^ l^^f^^^'^idTV + 9TAl5)b\ B{p)) . (32) 

The necessary matrix elements are given in Eq. (pT|l with Ti(0) = T2(0) for real photon 
emission. The calculation of decay rate gives, 

TMB - K*j) = {gly + gl^) (j^^^) <(1 - k*^fTl{Q). (33) 

The process B K*7 has been observed with a branching ratio 

Bexp{B K*-f) = (3.92 ± 0.20 ± 0.24) x 10"^ (34) 

Under the assumption Tnp{B K*'j) = Texp{B -f^*7), we get 

9tv + 9ta = ^-^^-oil X 10^'. (35) 

Comparing this constraint with those in Eq. we see that the process B — > K*'y puts 
a much stronger constraint on the g^y in comparison to that from B —>■ {K, K*)l^l~ . We 
substitute the above limit in Eq. (fT^ along with the approximation Ptv — Pta = 0.33 
GeV~^. The phase space integrals Itv and It a are essentially equal to each other. For 
electrons their value is 64 and for muons their value is 22. Then we get the branching ratios 
to be 

Bnp{B, ^ e+e-7) = 2.7ltlf, X 10-^ 

BNPiB, ^ /i+/i-7) = 9.18tlit X 10^'°. (36) 

These values are of the same order as SM predictions. Thus the stronger constraint on 
tensor/pseudo-tensor couplings coming from the experimentally measured value of B{B ^ 
K*'y) doesn't allow an enhancement of Bj^p{Bs 

Conclusions . The quark level interaction b s/+/~ is responsible for the three types 
of decays (a) semi-leptonic B — > {K, K*)l~^l~ , (b) purely leptonic Bg and also (c) 

leptonic radiative Bg — > I'^l^'y. It was shown in previously jlj| that if the purely leptonic 
branching ratio B{Bs Ai^Ai^)> 10^^ then new physics operators responsible for this have 
to be of the form scalar/pseudoscalar. Here we have shown that such scalar/pseudoscalar 
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operators have no effect on tlie leptonic radiative modes Bs — > Regarding otlier types 

of new physics operators, the vector/axial- vector operators can not enhance the branching 
ratios of Bg much beyond their SM values, given the constraints coming from 

the measured semi-leptonic rates. New physics operators in the form of tensor/pseudo- 
tensor also can not enhance Bg — > I'^l'j branching ratios given the constraints coming from 
B K*j. Thus we are led to the conclusion that the present data on 6 — > s 
transitions allow a large boost in B{Bs — > but not in B{Bs — > 

Acknowledgement 

This work grew out of discussions during the Ninth Workshop on High Energy Physics 
Phenomenology (WHEPP9) at Bhubaneshwar, India. We thank the members of the B- 
physics working group, especially Prof. David Hitlin, California Institute of Technology, 
USA, for valuable suggestions. We also thank Dr. Naveen Gaur, Delhi University, for an 
enlightening comment. 



[1] A. K. Alok, S.Uma Sankar, Phys. Lett. B 620 (2005) 61. 

[2] BaBar Collaboration: B. Aubert et ai, Phys. Rev. Lett. 91, (2003) 221802. 

[3] Belle Collaboration: A. Ishikawa et ai, Phys. Rev. Lett. 91, (2003) 261601. 

[4] G. Eilam, C. D. Lu, D. X. Zhang, Phys. Lett. B 391 (1997) 461. 

[5] T.M. Aliev, A.Ozpineci, M.Savci, Phys. Rev. D 55 (1997) 7059. 

[6] C.Q. Ceng, C.C. Lih, Wei-Min Zhang, Phys. Rev. D 62 (2000) 074017. 

[7] Y. Dincer, L.M. Sehgal, Phys. Lett. B 521 (2001) 7. 

[8] F. Kruger, D.Melikhov, Phys. Rev. D 67 (2003) 034002. 

[9] S. De Scotes-Genon, C. T. Sachrajda, Phys. Lett. B 577, (2003) 213. 
[10] D. MelikhoY, N. Nikitin, Phys. Rev. D 70 (2004) 114028. 
[11] CP. Korchemsky, D. Pirjol, T.M. Yan, Phys. Rev. D 61 (2000) 114510. 
[12] N. Gaur, S. Rai Choudhury, Phys. Lett. B 451, 86 (1999). 

[13] R. Forty, Pramana 63 (2004) 1135; N. Harnew, Nucl. Phys. Proc. Suppl. 120 (2003) 305. 

[14] W. Jaus, D. Wyler, Phys. Rev. D 41 (1990) 3405. 

[15] N.G. Deslipande, J. Trampetic, Phys. Rev. Lett. 60 (1988) 2583. 

[16] A. Ah, P.Ball, L.T.Handoko, G.Hiller, Phys. Rev. D 61, (2000) 074024. 



11 



[17] L. Lellouch, |hep-ph/0211359l D. Becirevic, |hep-ph/02113"40 



[18] BaBar Collaboration: B.Aubert et al, Pliys. Rev. D 70, (2004) 112006. 



